The Chaplygin dark star 
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We study the general properties of a spherically symmetric body described through the generalized 
Chaplygin equation of state. We conclude that such object, dubbed generalized Chaplygin dark star, 
should exist within the context of the generalized Chaplygin gas model of unification of dark energy 
and dark matter, and derive expressions for its size and expansion velocity. A criteria for the survival 
of the perturbations in the GCG background that give origin to the dark star are developed, and 
its main features are analyzed. 
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I. INTRODUCTION 
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The generalized Chaplygin gas (GCG) model has lately 
drawn some attention, mainly because it allows for a uni- 
fied description of dark matter and dark energy, the first 
dominating at early times and gradually transferring en- 
ergy to the dark energy component [1, 2]. The GCG 
model is consistent with various classes of cosmological 
tests, such as the Cosmic Microwave Background Ra- 
diation [3], supernovae [4], gravitational lensing [5] and 
gamma-ray bursts [6]. As with other competing candi- 
dates to explain the overwhelming energy density of the 
present Universe, GCG is naturally constrained through 
cosmological observables. 

It is quite interesting that the GCG equation of state 
is that of a polytropic gas [9] , although one with a neg- 
ative polytropic index. This hints that one could look 
for astrophysical implications of the model, and hence 
hope for yet another approach to the problem of con- 
straining the allowed space for its parameters (see, e.g. 
Ref. [10]). In this work we argue that a GCG dark star 
may arise from a density fluctuation in the cosmological 
GCG background. In what follows we shall character- 
ize these objects, look at their evolution and account for 
their initial probability of appearance within the GCG 
background. 
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where a is the scale factor of the Universe and B a pos- 
itive integration constant. This result shows a strik- 
ing property of the GCG, namely that at early times 
it behaves as non-relativistic dark matter {p c h °c a -3 ), 
while at late times it acts as a cosmological constant 
{pch — const.). One can algebraically decompose the 
GCG into a dark matter and a dark energy component, 
which evolve such that the transferring of energy occurs 
from the former to latter [7]. This can be used to show 
that, while the dark energy component is spatially homo- 
geneous, the dark matter component allows for structure 
formation at early times, when it dominates [2, 7, 8]. 

For convenience, one defines the parameter A s = 
A/p^uq, where Video is the dark energy density, p cr o the 
critical density and p c ho the GCG energy density, all at 
the present. Assuming, as observations suggest, the con- 
dition fldmo + ^rfeo = lj where fldmo is the dark matter 
density (dropping the small baryon contribution), and 
taking ciq = 1, yields the constraint B = QdmaPcrOp^hO- 



III. POLYTROPIC STARS 



II. THE GENERALIZED CHAPLYGIN GAS 



The GCG model is based on the equation of state 



In order to deal with stellar structure in general rela- 
tivity, one considers that the spherical body behaves as 
a perfect fluid, characterized by the energy-momentum 
tensor 
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where A and a are positive constants and < a < 1 
(see however Ref. [11] for reasons to consider a > 1); 
the negative pressure hints that the GCG is related to a 
cosmological constant. The case a = 1 corresponds to the 
Chaplygin gas [12]. In a Friedmann-Robertson- Walker 
cosmology, the relativistic energy conservation yields 



= {p + Pju^u" + Pg^ 



(3) 



where g^ v is a spherically symmetric Birkhoff met- 
ric. The Bianchi identity implies that V^T^ = 0, 
from which follows the relativistic Tolman-Oppenheimer- 
Volkov equation [9] , 
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where m(r) = 4ir p(r')r' 2 dr' . This equation collapses 
to the classical Newtonian hydrostatic equilibrium equa- 
tion 
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if the following conditions are satisfied: 

Gm(r)/r < 1 , 
47rr 3 P(r) < m(r) , 
P(r) < p(r) . 



(5) 



(6) 
(7) 
(8) 



The polytropic gas model for stellar like structure as- 
sumes an equation of state of the form P = Kp n+1 / n , 
where n is the polytropic index, which defines interme- 
diate cases between isothermic and adiabatic thcrmody- 
namical processes, and K is the polytropic constant, de- 
fined as 



K = N n GM^ n -^/ n R^- n ^ n 
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with 
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where R is the star's radius, M its mass and £i, defined 
by = 0, corresponds to the surface of the star (cf. 

below) . Actually, this definition states that all quantities 
tend to zero as one approaches the surface. 

This assumption leads to several scaling laws for the 
relevant thermodynamical quantities, 
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where p c , T c and P c are the density, temperature and 
pressure at the center of the star [9]. Notice that the 
scaling law for temperature requires the assumption that 
the gas behaves as an ideal one. This is not the case of 
the GCG. 

The function 9 is a dimensionless function of the di- 
mensionless variable £, related to the physical distance 
to the star's center by r = /?£, where 
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The function 6(£) obeys a differential equation aris- 
ing from the equilibrium condition of Eq.(5), the Lane- 
Emden equation: 
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Notice that the physical radius and mass of the spherical 
body appear only in the polytropic constant, and hence 
the behavior of the scaling function is unaffected by 
these. Therefore, the stability of a star is independent of 
its size or mass, and different types of stars correspond 
to different polytropic indices n. This scale-independence 
manifests in the symmetry of the Lane-Emdcn equation 
(15) under homology transformations. The first solar 
model ever considered, developed by Eddington in 1926, 
was that of an n — 3 polytropic star. Although somewhat 
incomplete, this simplified model gives rise to relevant 
constraints on the physical quantities. 

In what follows, we shall use the Lane-Emden equa- 
tion to derive the properties of a generalized Chaplygin 
dark star, given that conditions (6)-(8) are shown to be 
fulfilled. This enables the use of the Newtonian approxi- 
mation (5), which asides its simplicity allows for a prompt 
interpretation of the GCG as a polytropic gas subject to 
the Lane-Emdcn equation of motion. The generality of 
this procedure can be used in various cases of physical in- 
terest, as for instance, when studying the effect of scalar 
fields on the stellar structure [10]. 



IV. THE GENERALIZED CHAPLYGIN DARK 
STAR 

As already discussed, the GCG model is cosmological 
in nature, and most bounds on the parameters a and A s 
are derived from cosmological tests. However, a quick 
look at the GCG equation of state (1) indicates that the 
it corresponds to a polytrope with a negative polytropic 
constant and a negative pressure. At first glance, this 
seems to indicate that no valid analysis of a GCG at 
an astrophysical context can proceed, since a spherical 
body constituted by such exotic gas would experience an 
outward pressure that would prevent it from being sta- 
ble. However, the following argument shows that such 
an objection is circumvented by the presence of the cos- 
mological GCG background. 

The first logical step for the construction of a symmet- 
ric body with the GCG equation of state should be, as 
for all polytropes, the solution of the related Lane-Emden 
equation. Firstly one notes that this stems from the hy- 
drostatic equation (5). As already seen, this is directly 
derived from the general relativity equations, assuming 
the energy-momentum tensor of a perfect fluid: no as- 
sumption whatsoever is made concerning the pressure or 
density, nor the equation of state relating these quan- 
tities. Hence, one can use this equation and, through 
the usual derivation, the related Lane-Emden equation; 
as stated before the only concern is if one can neglect 
the higher-order relativistic terms present in Eq.(4), thus 
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working in the Newtonian limit. This will be explicitly 
shown in the case under investigation. 

The polytropic equation of state P — p n+1 / n shows 
that the GCG can be assigned a negative polytropic in- 
dex n = — 1/(1 + a). Next, a comparison with Eq. (1) 
yields K = —A, since the pressure is negative. This re- 
quires some caution: indeed, the direct application of 
the coordinate transformation between the physical ra- 
dial coordinate r and the coordinate £ given by by r = 
with j3 defined in Eq. (14) yields 
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where D = [pi 1 n)/n /AttG} 1 / 2 ; since a > and K < 0, 
the above quantity is imaginary. To avoid this, we define 
the coordinate £ through the same equation r = /3£, but 
with K replaced by \K\ = A > in Eq. (14), obtaining 
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The negative sign of K will, of course, manifest itself 
in the terms of the Lane-Emden equation; explicitly, one 
gets 
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As in the usual Lane-Emden equation, one has as 
boundary conditions 0(0) = 1 and 6"(0) = 0. This gives 
rise to a positive derivative for £ > 0, indicating that 9 
is a smoothly increasing function. This is related not to 
the negative polytropic index, but to the negative pres- 
sure of the GCG; as a consequence, Eqs. (11) and (13) 
indicate that the pressure inside the dark star increases 
(in absolute value), while the density decreases. This is 
key to our study, since it shows that a GCG spherical 
body accretes, as expected for a star. 

In the usual Lane-Emden equation, the criteria con- 
cerning the size of a star is given by #(£i) = 0, corre- 
sponding to a surface of zero density, pressure and tem- 
perature. In a GCG dark star, the question is more con- 
voluted: a vanishing density yields infinite pressure (and 
conversely), which are rather unphysical choices for the 
boundary of any astrophysical object. Furthermore, since 
the function 9 is increasing, the density p oc 6 n vanishes 
at an infinite distance, while the pressure P oc —9 1+n 
does not vanish at all. 

As a solution for this issue, one recalls that the GCG 
object is embedded in a cosmological background. Hence, 
a GCG dark star should not be taken as an isolated body, 
but rather as a spike on the overall cosmological back- 
ground of density p c h and negative pressure P c h- There- 
fore, its boundary should be signalled by the matching of 
the inner and outer pressures and densities, as indicated 
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FIG. 1: The function 0(£)/8 for a relative density 8 = 10, 
100, 1000 (dashed, dot-dashed and full lines respectively) as 
a function of £/5, assuming a = 0.2 and A s = 0.7. 



in Fig. 1. Both conditions are, of course, equivalent, 
given the GCG equation of state (1). From Eqs. (11) 
and (13), this equates to 



6>(6) = <T V " = S 1+c 
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where one defines the ratio between central and the back- 
ground density, S = p c / Pch- Hence, one gets a correspon- 
dence between the central density of the dark star, that 
is, the height of the energy density fluctuation, and its ra- 
dius. This argument shows that the Chaplygin dark star 
cannot be taken merely as an isolated body having a com- 
mon equation of state with the GCG, but must be viewed 
instead as a perturbation to the flat GCG background. 
This is advantageous, since one can use the constraints 
available on the GCG model to ascertain its properties. 
The only constraint affecting this quantity is S » 1, since 
we assume that the density perturbation must be large 
enough so we can view it as a physical object, not merely 
a small fluctuation on the GCG background. 

Notice that, since the scaling function is com- 
pletely specified by the two boundary conditions 9(0) = 1 
and 9'(0) = 0, neither p c nor p c h affect each other: p c is 
"scaled out" of the problem through the definition (11) 
and(13) and p c t merely sets the criteria for the surface 
of the star, through Eq. (19). Hence, although p c h varies 
with time, there is no contradiction in assuming a con- 
stant central density p c ; this simplifies our study, since 
Eq. (17) then yields a constant coordinate scaling coeffi- 
cient (3. Furthermore, given that the cosmological back- 
ground density p c h is decreasing, a constant central den- 
sity p c indicates that the density ratio 5 increases and, 
therefore, £i expands towards a final radius = /3£ioo, 
acquiring mass in the process. One can argue that, due 
to energy conservation, the background density should 
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decrease in order to compensate this, but this is a minor 
effect that can be neglected. Also, since the Chaplygin 
dark star is not an isolated object, it is reasonable to 
assume that neither its mass nor radius should be held 
constant, but instead must vary as it dilutes itself on 
the overall cosmological background; instead, the central 
density p c arises as the natural candidate for distinguish- 
ing between these objects. 

The absolute magnitude of the perturbation results 
from the dynamics ruling the generation of a perturba- 
tion, via a probability law arising from the fundamen- 
tal physics underneath the GCG model; this, of course, 
should naturally disfavor very large perturbations. Fur- 
thermore, since any relative perturbation to the homoge- 
neous energy density profile is local, its occurrence should 
not depend on cosmological quantities and the probabil- 
ity distribution should depend only on the relative per- 
turbation 5, and not explicitly on the scale factor a. A 
putative candidate could be a normal probability distri- 
bution given by 



f(5) = f exp[-gS 2 



(20) 



where /o is a normalization factor, g a parameter de- 
pendent on A s and a. Of course, more complicated ex- 
pressions for f(5) are possible, depending on the inner 
workings of the fundamental physics behind the GCG 
model. 

The expansion velocity of a dark star can be shown to 
be given by 
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where the prime denotes differentiation with respect to £ 
and H = a/a is the rate of expansion of the cosmological 
background. The Friedmann equation allows one to write 
the latter as 
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One can see that, as p c h approaches a constant value at 
late times, the expansion tends to zero. Also, one can 
derive the dependence of v e on 8 by noticing that the 
scaling coefficient (3 runs with 5~( 1+a / 2 ^ and the term in 
brackets in Eq. (21) scales with 5, while is found 

numerically to always be of order unity. Hence, one 
concludes that the expansion velocity depends weakly 
on S, v e cx <5~ Q / 2 ~ <5~ al , given the chosen value of 
a = 0.2 [7, 13]. By the same token, since numeri- 
cally one finds that £i cx S, it can be concluded that 
n = /3£i oc 5- a ' 2 ~ 5~ 01 . 

Given that the GCG tends to a smooth distribution 
over space, most density perturbations tend to be flat- 
tened within a timescale related to their initial size and 
the characteristic speed of sound v s = (dP/dp) 1 / 2 . Inside 



the dark star, the equation of state (1) and the definitions 
(11) and (13) yield 



aAm_ , 0(0 (Pern 
which at the surface amounts to 
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One sees that the maximum sound velocity occurs at the 
surface of the star; since a < 0.6 and 0.6 < A s < 0.8 (see 
first references in [3]), this should be smaller than the 
present value of v s ^ max — 0.693 (in units of c). A plau- 
sible criteria for the survival of an initial perturbation is 
given by v SiC h < v e ,o, where w Ci o is the initial expansion 
velocity. Equating Eqs. (21) and (24) yields, after a little 
algebra 
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Numerically, one finds that the left-hand side of Eq. 
(25) is approximately constant. At early times, the 
GCG behaves as cold dark matter, with p c h oc a~ 3 and 
H oc a~ 3 / 2 , and the right hand side of condition (25) is 
constant. At late times, when the GCG acts as a cosmo- 
logical constant, both p c h and the expansion rate H are 
constant, and the rhs then scales as a~ 3 ( 1+a \ and thus 
decreases with cosmic time. Thus, most dark stars are 
created at early times; this is consistent with the usual 
interpretation of the GCG as a model where dark matter 
dominates at early times, allowing for structure forma- 
tion, while at late times the dark energy component takes 
over. 

Eq.(24) and the GCG equation of state (1) allows one 
to rewrite condition (8) as 
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where we have used p(r) > p c h > Pchoo for any redshift 
z, with pchoo = A 1 / 1+Q the limit for the background cos- 
mological density when a — > oo. Relativistic corrections 
could be important if the pressure and density are of 
the same order of magnitude. However, one can use the 
bound 
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to ascertain that, for redshifts typical for structure for- 
mation, z = z c = 15 and a set of GCG parameters 
a = 0.2, A s = 0.7, one gets \P(r)/p(r)\ < 10~ 4 ; higher 
redshifts provide an even lower upper bound. For a 
much recent z = 1, the same set of parameters yields 
\P(r)/p(r)\ < 0.188, which still validates the Newtonian 
approximation, although to a lesser extent. This is not 
troublesome, since most dark stars are assumed to nucle- 
ate at an early age. Also, since the above condition only 
provides an upper limit for P(r)/p(r), a more complete 
calculation can still validate the Newtonian approxima- 
tion, depending on the value of 5. 

Assuming that all dark stars have expanded up to their 
final size (which follows from the stabilization of the GCG 
as a dark energy), one can write the mass contribution 
of those created when the Universe had a size a(t): 

-^=J o J o p(r,Pc)f(Pc)r 2 drdp c (28) 

fPn ril(Pa) 

=/? 3 / / PcO n {^pc)f{Pc)e^dp c , 

Jo Jo 

where n = and the dependence on the integration 
variables is made explicit. Integrating over time one gets 
the mass contribution of all generations of dark stars, 
M DS = J a ° M'(a)da, where M'(a) = dM/da. 

A comparison between known observational bounds 
and numerical integration of the above results can then 
be used to constraint the GCG parameters A s and a, 
namely through supernovae data, gravitational lensing 
results and other dark matter searches. This will be con- 
sidered elsewhere. 



V. NUMERICAL RESULTS 

In order to substantiate our arguments, in this section 
we present some numerical examples; we shall study the 
proposed scenario for the "typical" values of the GCG 
model; one takes a = 0.2 and A s = 0.7. A future study 
based on this results could embrace a wider range of pa- 
rameters and probe the creation of dark stars at early 
stages of the Universe, providing further refinement to 
the already known bounds (see, e.g. Ref. [13] for a sum- 
mary of the existing constraints). 

A numerical integration of the modified Lane-Emden 
equation (18) produces the results plotted in Fig. 1. In 
Table I we draw different scenarios, in order to ascer- 
tain the dimensions of dark stars nucleated at different 
ages of the Universe. One can see that at a redshift of 
z = z c = 15, presumably typical for structure formation, 
even a small perturbation 6 = 5 produces an overwhelm- 
ingly large object, with about 3000 times the mass and 
20 times the diameter of the Milky Way. Since these 
dimensions scale with (3, which decreases with p c , one 
probes higher redshifts in order to obtain smaller dark 
stars. Therefore, at z = 50 and S = 10, one obtains 



an object with approximately the size and double the 
mass of our galaxy. A larger perturbation S = 100 yields 
approximately the same size, but a ten-fold increase in 
mass. 

Going further back in time, a darkstar born at z = 100 
with S — 10 (100) has about one-hundredth (one-tenth) 
the mass of the Milky Way and one-tenth its diame- 
ter. Finally, 6 = 100 and an extremely high redshift 
of z = 500, deep within the so-called "dark ages", yield 
a dark star with 1.6 x 10 6 solar masses and a radius of 
7.8 pc, dimensions similar to those ascribed to super- 
massive black holes in active galactic nuclei. 

The above discussion is by no means definitive, and 
only serves to illustrate the concept developed in this 
study. Nevertheless, one might be surprised by the un- 
physically large size of a dark star hypothetically nu- 
cleated at the redshift typical for structure formation, 
z c = 15. However, notice that this describes the con- 
densation of baryonic matter interacting gravitationally 
with dark matter. The dark star scenario poses quite a 
different mechanism, where the GCG (in an era where 
its dark matter component dominates) is the sole con- 
stituent of the spherical body. Hence, it is reasonable to 
assume that bodies of astrophysical dimensions can arise 
much earlier in the history of the Universe. A precise 
description would of course imply the nucleation prob- 
ability distribution f(8), since this is the fundamental 
quantity ruling the onset of perturbations nucleation on 
the GCG background. 

To ascertain the validity of the Newtonian limit from 
which the Lane-Emden equation is derived, conditions 
(6) and (7) can now be checked with a simple calcula- 
tion. An inspection of the Table I shows that £i is of the 
order of 5. Also, Fig. 1 shows that while slowly 
increases and is of order unity, the scaling function #(£) 
grows regularly; hence, one can use the power-law ap- 
proximation ~ at; b , where a and b are coefficients of 
order unity, since it does not introduce large deviations 
from a full numerical calculation. Using Eq. (11), this 
yields 



mW ~ JTnb p{r)r3 ' (29) 
and condition (6) becomes 



— b Gp(ry « 1 . (30) 

Since p(r) oc r nh , the left-hand side of Eq. (30) scales 
with r 2+nb and, since \nb\ ~ 1, is an increasing func- 
tion. Therefore, it is majored by its value at r = n, 
amounting to ~ AnGpchrf. Using Table I one finds that 
it attains a maximum value of 5.12 x 10~ 5 for S = 5, 
z = z c = 15, thus concluding that condition (6) is veri- 
fied. Notice that, for a fixed redshift, this ratio is approx- 
imately independent of 5. This is due to the very weak 
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5 = 5, z = 15 




8 = 10 , z = 50 




8 = 100 , 2 = 50 




8 = 10 , z = 100 




8 = 100 , z = 100 




8 = 100 , z = 500 




9.04 




19.1 




239 




19.1 




239 




239 


#'(6) 


0.850 




0.916 




1.14 




0.916 




1.14 




1.14 


p c (Kg.rn~ 3 ) 


6.11 x 10" 23 




3.96 x 10~ 21 




3.96 x 10~ 20 




3.07 x 10" 20 




3.07 x 10~ 19 




3.75 x 10~ 17 


P (pc) 


7.59 x 10 4 




772 




61.4 




81.0 




6.44 




3.26 x 10~ 2 


ri (pc) 


6.86 x 10 5 




1.48 x 10 4 




1.47 x 10 4 




1.55 x 10 3 




1.54 x 10 3 




7.81 


M/M Q 


1.73 x 10 15 




1.13 x 10 12 




1.11 x 10 13 




1.02 x 10 10 




9.97 x 10 10 




1.58 x 10 6 




5.12 x 10~ 5 




7.66 x 10" 7 




7.58 x 10~ 7 




6.55 x 10" 8 




6.48 x 10~ 8 




2.03 x 10 -10 




1.89 x 10~ 2 




2.33 x 10" 3 




2.35 x 10" 3 




6.81 x 10~ 4 




6.86 x 10~ 4 




3.84 x 10" 6 




5.09 x 10~ 3 




6.32 x 10~ 4 




6.32 x 10~ 4 




1.85 x 10~ 4 




1.85 x 10~ 4 




1.03 x 10" s 


6>'(£i)<5~ a/2 


0.724 




0.728 




0.722 




0.728 




0.722 




0.722 



TABLE I: Numerical results for the quantities £i, 0'(£i), p c , P, ri, M, 4TrGp ch rf /c 2 , v e0 , v s and 0'(£i)<5~ a/2 , for 8 = 10, 100 
and redshifts z = 15, 50, 100, 500. 



scaling of the physical radius with the relative perturba- 
tion, n oc 5~ a / 2 ~ 5~ 0A , for the chosen value a — 0.2. 

In a similar fashion, the power-law approximation al- 
lows one to write condition (7) as 

P(r)«P(r)^ b , (31) 

which, since \nb\ ~ 1, is equivalent to condition (8) and 
hence also satisfied at early ages. Thus, the Newtonian 



approximation implicit in the Lane-Emden equation is 
valid for the cases studied and, given the smallness of 
the values encountered in its evaluation, it is also appli- 
cable in a broader range of the nucleation redshift z and 
relative densities 5. 

Given the indicated values for the initial expansion ve- 
locity v e o and the surface sound velocities v s (which de- 
pends only on the redshift) the inequality (25) is valid 
for the chosen values, and thus the corresponding dark 
stars do not collapse at birth. 



VI. CONCLUSIONS 

In this study we have analyzed the properties of spher- 
ical bodies with a polytropic equation of state of nega- 
tive index, in the context of the GCG dark energy/dark 
matter unification model. We have considered the associ- 
ated Lane-Emden equation and looked at the qualitative 
behavior of its solution; amongst the results we find the 
conditions for fluctuations to be attenuated or to develop 
as dark stars. Our criteria is based on the condition that 
the sound velocity does not exceed the expansion veloc- 
ity of the dark star when it nucleates v s , c h < v e ,o- This 
enables the computation of the mass contribution of the 



dark stars at present times, which can then be used to 
constraint the GCG parameters A s and a, providing an- 
other testing ground for this fascinating model. 

Note added: While finalizing this work we became 
aware of the study of stable dark energy objects [14] 
and halos of fc-essence [15]. Even though the motiva- 
tion of both works are somewhat similar, our approaches 
are quite different. 
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